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Nuclear composition and heating 
in accreting neutron-star crusts 
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Abstract. Nuclear reactions in accreting neutron-star crusts and the heat release accompanying them are studied, 
under different assumptions concerning the composition of the outermost layer formed of the ashes of X-ray bursts. 
Particular examples of ashes containing nuclides with A ~ 90 — 110 are considered and compared with a standard 
j4 = 56 case. In all cases, evolution of a crust shell is followed from 10* g cm~^ to a few times 10^^ g cm~^. The 
total crustal heating produced in the non-equilibrium processes in the accreting crust is 1.1 — 1.5 MeV per one 
accreted nucleon. The composition of the accreted crust at densities exceeding the threshold for the pycnonuclear 
fusion (p > 10^^ g cm~^) is essentially independent of the assumed initial composition of the X-ray burst ashes. 
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1. Introduction 

Accretion of matter onto neutron stars in close binary sys- 
tems plays crucial role in many objects studied by the X- 
ray astronomy, such as X-ray bursters. X-ray pulsars, and 
transient X-ray sources. Particular attention was focused 
recently on soft X-ray transients (SXRTs) in quiescence, 
when the accretion from a disk formed of matter flow- 
ing from the low-mass companion star is switched off or 
strongly suppressed. The nature of the X-ray radiation 
during quiescence is still a matter of debate. Typically, 
the quiescent emission is much higher than the expected 
one for an old cooling neutron star. It has been suggested 
that this is due to the fact that the interiors of neutron 
stars in SXRTs are heated-up, during relatively short pe- 
riods of accretion and bursting, by the non-equilibrium 
processes associated with nuclear reactions taking place in 
the deep layers of the crust (Brown et al. 1998). This addi- 
tional crustal heating, combined with appropriate models 
of neutron-star atmosphere and interior, can be used to ex- 
plain observations of SXRTs in quiescence. Dependence of 
the luminosity in quiescence on the structure of neutron- 
star core, and particularly on the rate of neutrino cooling, 
opened a new possibility of exploring the internal struc- 
ture and equation of state of neutron stars via confronta- 
tion of theoretical models with observations of quiescent 
SQRTs (see, e.g., Ushomirsky & Rutledge 2001, Colpi et 
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al. 2001, Rutledge et al. 2002, Brown et al. 2002, and ref- 
erences therein). 

The crustal heating used in the SQRTs models was 
described using the results of the calculations of Haensel 
& Zdunik (1990) (hereafter referred to as HZ). In the HZ 
model, an outer layer of neutron star, formed in X-ray 
bursts, is composed of ^^Fe. Such a composition, corre- 
sponding to the maximum nuclear binding, was predicted 
by the standard models of the thermonuclear X-ray bursts 
of that time. Then, the ^^Fe layer is sinking in the neutron 
star interior under the weight of accreted matter. Under 
increasing pressure, the composition of matter is chang- 
ing in a sequence of nuclear reactions: electron captures, 
neutron emission and absorption, and finally, at densities 
exceeding 10^^ g cm^"^, also pycnonuclear fusion. As the 
nuclear processes are proceeding off-equilibrium, they are 
accompanied by the heat deposition in the crustal matter. 
For the HZ model, the total crustal heating amounts to 
1.4 MeV per one accreted nucleon. This heating is mostly 
supplied by the pycnonuclear fusion processes in the inner 
crust at p = 10^^ — 10^^ g cm^'^. 

Very recently, new X-ray burst simulations, with a net- 
work of nuclei much larger than the previously used net- 
works which ended with the iron-peak nuclei, were carried 
out by Schatz et al. (2001). Extending the network of pos- 
sible nuclear reactions resulted in nuclear ashes composed 
predominantly of nuclei much heavier than iron, with A 
up to ~ 112. The question arises how changing the ini- 
tial composition from ^^Fe to, e.g., ^"''Te will change the 
crustal heating and the composition of the deeper layers of 
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the inner crust. The present Letter is devoted to answering 
this question. 

Brief description of the dense-matter model and of the 
processes taking place in an accreting neutron-star crust is 
given in Sect. 2. Results for the composition and heating, 
obtained for various initial compositions and dense-matter 
models, are presented in Sect. 3. Final Sect. 4 is devoted 
to a discussion of our results and conclusion. 

2. Non-equilibrium nuclear processes 

In what follows we briefly describe the nuclear evolution 
scenario of HZ as applied in the present study (earlier 
studies on various aspects of the non-equilibrium neutron- 
star crusts were done by Vartanyan & Ovakimova 1976, 
Bisnovatyi-Kogan & Chechetkin 1978, Sato 1979). Under 
conditions prevailing in accreting neutron-star crust at 
p > 10^ g cm~'^ matter is strongly degenerate, and is 
relatively cold (T < 10^ K), so that thermonuclear pro- 
cesses involving charged particles are blocked. At the den- 
sities lower than the threshold for the pycnonuclear fusion 
Ppyc ~ 10^^ g cm~^, the number of nuclei in an element of 
matter does not change during the compression resulting 
from the increasing weight of accreted matter. For sim- 
plicity, we assume that only one nuclear species {A, Z) is 
present at each pressure. Due to the nucleon pairing, sta- 
ble nuclei in dense matter have even N — A — Z and Z 
(even-even nuclides). In the outer crust, in which free neu- 
trons are absent, the electron captures which proceed in 
two steps, 

{A, Z) + e- {A, Z-l) + v, , 
[A, Z-l) + e- [A,Z-i) + j/e + Qc . (1) 

lead to a systematic decrease of Z (and increase of = 
A — Z) with increasing density. The first capture in Eq. ^ 
proceeds in a quasi-equilibrium manner, with a negligible 
energy release. It produces an odd-odd nucleus, which is 
strongly unstable in dense medium, and captures a sec- 
ond electron in an non-equilibrium manner, with energy 
release Qc- A fraction of this energy, Q^, is deposited in 
the matter and heats it, but most of it, Q^, is carried away 
by the neutrino, so that Qc = Qd + Qv The effective heat 
deposited in matter is 

Qd ^ ??(Me - A) + Qoxc , (2) 

where /Ug is the electron Fermi energy (including the rest 
energy) and A is the energy threshold for the first quasi- 
equilibrium electron capture. Qcxc is the excitation energy 
of the final nucleus , and t] ranging from 1/6 to 1/4 ac- 
counts for the neutrino energy losses. We put Qoxc — 
in the actual simulations, except for one case with ini- 
tial ^^Fe composition when Qcxc is experimentally known 
(see HZ). We have 77 1/6 for ^ A and 77 = 1/4 for 
(/Ue — A)//Ze <C 1. In the crust evolution simulation the 
latter condition is fulfilled and therefore we put r\ = 1/4. 



Above the neutron-drip point (p > pnd), electron cap- 
tures trigger neutron emissions, 

(A, Z) + e- — > (A, Z - 1) + , 

(A, Z - 1) -I- — > (yl - k, Z - 2) -I- k n z^e + Qc , (3) 

where the number of emitted neutrons "k" is even. Due 
to the electron captures, the value of Z decreases with in- 
creasing density. In consequence, the Coulomb barrier pro- 
hibiting the nucleus- nucleus reaction lowers. This effect, 
combined with the decrease of the separation between the 
neighboring nuclei, and a simultaneous increase of energy 
of the quantum zero-point vibrations around the nuclear 
lattice sites, opens a possibility of the pycnonuclear re- 
actions. The pycnonuclear fusion timescale Tpyc is a very 
sensitive function of Z. The chain of the reactions leads 
to an abrupt decrease of Tpyc typically by 7 to 10 orders 
of magnitude. Pycnonuclear fusion switches-on as soon as 
Tpyc is smaller than the time of the travel of a piece of mat- 
ter (due to the accretion) through the considered shell of 
mass A'fshcn(A^, ^), ^acc = ^'h\ic\\l^i- The masses of the 
shells are of the order of 10~^ Mq. As a result the point 
where the pycnonuclear reaction takes place is very well 
defined and the chain of reactions (jJl in several cases is 
followed by the pycnonuclear reaction on a timescale much 
shorter than Tacc- Denoting Z' = Z — 2, we have then 

(A, Z') + (A, Z') (2A, 2Z') + Qi , 
(2A, 2Z') (2A - k', 2Z') + k' n + Q2 , 

+ Q3 , (4) 

where ". . ." correspond to some not specified chain of the 
electron captures accompanied by neutron emissions. The 
total heat deposition in matter, resulting from a chain 
of reactions involving a pycnonuclear fusion, is Qpyc 
Qi + Q2 + Qs.d where label "d" indicates that only the 
fraction deposited in the matter is included. 

Our model of atomic nuclei is described in HZ. Using 
our nuclear-evolution code, we evolved an element of mat- 
ter consisting initially of nuclei (Ai,Zi) immersed in an 
electron gas, from pi = 10^ g cm~^ to p = pf > 5 x 
2^q13 g Q^j, results for the composition and crustal 

heating are presented in the next section. 

3. Composition and heating 

The compositions of accreted neutron star crusts are 
shown in Fig. and in Table ^ These results describe 
crusts built of accreted and processed matter up to the 
density 5 x 10^^ g cm~^. At a constant accretion rate 
M = Af_9 X 10^9 Mo/yr this will take - 10^ yr/M-g. 
During this time, a shell of X-ray burst ashes will be com- 
pressed from ^ lO'* g cm"'^ to ^ 10^^ g cm~'^. Two differ- 
ent compositions of X-ray bursts ashes at ^ lO'* g cm^'^ 
were assumed. In the first case, A\ — 56, and the HZ 
scenario is reproduced. In the second case, we consider 
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Fig. 1. Z and N of nuclei versus matter density in an 
accreting neutron-star crust. Solid line: A\ = 106; dotted 
line: A\ = 56. Each change of N and Z, which takes place 
at a constant pressure, is accompanied by a .jump in den- 
sity (see HZ for detailed discussion of this point). Small 
steep segments connect the top and the bottom density of 
thin reaction shell. Arrows indicate positions of the neu- 
tron drip point. 



an example of nuclear ashes obtained by Schatz et al. 
(2001). To be specific, we assume Ai — 106. The value 
of Zi — 46 stems then from the condition of beta equilib- 
rium at p = 10^ g cm~^. The compositions in the outer 
crust, where the only processes are the electron captures, 
are strongly influenced by the initial conditions. Up to 
the neutron-drip point, the difference by a factor of about 
two between the values of Z and N for Ai = 106 and 
Ai — 56 is conserved. It should be noted that in the case 
Ai = 106 the number of beta captures in the outer crust 
is about 2 times larger, but each reaction is accompanied 
by the density jump ~ 5%, about half of these in the 
case Ai — 56 (we have similar situation for the energy 
release in a single shell). One notices a dramatic effect 
of the neutron drip triggered by an electron capture at 
p = PND- We get pND = 6 X 10^^ g cm~^ for Ai — 56 and 
Pnd = 8 X 10^^ g cm"'^ for Ai = 106. Just after pnd both 
N and Z of nuclei decrease in a long chain of the neutron 
emissions followed by the electron captures. After the py- 
cnonuclear fusion is switched-on at Ppyc — 10^^ g cm~^, 
the two compositions converge, and stay very close, up 
to the largest densities beyond which the validity of the 
HZ model becomes questionable. We checked that this is 
a generic property of the Z, N evolution, which does not 
depend on specific values of Ai, Zi, or on the details of the 
nuclear model used (see Sect.0J. 



Fig. 2. Heat per one accreted nucleon, deposited in the 
crust, for two models with different initial A. Solid vertical 
lines (ended with circles): Ai = 106; dotted lines (ended 
with crosses): Ai — 56. Vertical lines are positioned at the 
densities at the bottom of the reaction shell. 

In Fig. 121 we show the heat deposited in the matter, 
per one accreted nucleon, in the thin shells in which non- 
equilibrium nuclear processes are taking place. Actually, 
reactions proceed at a constant pressure, and there is a 
density jump within a thin "reaction shell" . The vertical 
lines whose height gives the heat deposited in matter are 
drawn at the density at the bottom of the reaction shell. 

One notices a specific dependence of the number of 
heat sources and the heating power of a single source 
on assumed Ai. Let us start with the outer crust (up- 
per panel of Fig. In the case of Ai = 56 the number 
of sources is smaller, and their heat-per-nucleon values 
q are larger, than for Ai = 106. Hovever, the total de- 
posited heat-per-nucleon is quite similar, 0.041 and 0.039 
MeV/nucleon for Ai — 56 and Ai = 106, respectively. 
Similar features are seen in the inner crust (lower panel 
of Fig. 121). The total crustal heating is 1.54 MeV/nucleon 
and 1.12 MeV/nucleon for Ai — 56 and Ai = 106, respec- 
tively. The difference ~ 0.4 MeV/nucleon between these 
two cases is mainly due to an additional pycnonuclear 
reaction (the first one) in the case Ai = 56, which re- 
sults in the convergence of two evolutionary scenarios at 
p ~ 2 X 10^^ g cm~'^. This pycnonuclear fusion is accom- 
panied by the larger energy release than the subsequent 
beta captures and neutron emissions in the case Ai — 106. 
The nearly exact convergence of the cases Ai — 56 and 
Ai = 106 for p > 10^^ g cm~'^ is connected with the fact 
that heavier nucleus has iVi and Zi which are nearly dou- 
ble of those of ^^Fe. In the case of the initial nuclei be- 
tween Ai = 56 and Ai = 106 the situation is similar with 
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Table 1. Non-equilibrium processes in the crust of an accreting neutron stars assuming that the X-ray ashes consist 
of ^"^^Pd. P and p are pressure and density at which the reaction takes place. Ap/p is relative density jump connected 
with reaction, q is the heat deposited in matter, X„ is the fraction of free neutrons among nucleons, in the layer just 
above the reaction surface. 



p 


p 


reactions 


X„ 


Ap/p 


q 


(dyn cm~^) 


(e cm"*) 






% 


(keV) 


9.235 X lO^'' 


3.517 X 10"* 


'""Pd Ru - 2e~ + 2ue 





4.4 


5.7 


3.603 X 10^'^ 


5.621 X 10°^ 


""'Ru ^'"'' Mo - 2e" -1- 2z/e 





4.6 


5.7 


2.372 X 10^* 


2.413 X 10'° 


106|^Q ^106 _ 2e~ + 2l>e 





4.9 


5.6 


8.581 X 10^^ 


6.639 X 10'° 


i06zr ^'°« Sr - 2e- + 2i/e 





5.1 


5.6 


2.283 X 10^^ 


1.455 X lO" 


106g^ ^106 j^j. _ 2g- + 2Ue 





5.4 


5.5 


5.025 X 10^^ 


2.774 X 10'' 


"'6Kr ^'°6 Se - 2e- + 2u, 





5.7 


5.5 


9.713 X 10^" 


4.811 X 10" 


"'«Sc^'°6 Ge-2e- +2ue 





6.1 


5.5 


1.703 X 10^° 


7.785 X 10" 


'"''Ge ->°^ Ni + Un - Ae' + 4i/e 


0.13 


13.2 


77.6 


1.748 X 10*° 


8.989 X 10" 


92^1 ^86 -pe + 6n- 2e~ + 2ve 


0.19 


6.9 


39.2 


1.924 X 10^° 


1.032 X 10'^ 


86pg ^80 Cj. g„ _ + 2Ve 


0.25 


7.3 


43.1 


2.135 X 10^° 


1.197 X 10'^ 


™Cr -^''^ Ti -h 6n - 26" + 2v^ 


0.30 


7.7 


47.4 


2.394 X 10*° 


1.403 X 10'^ 




0.36 


8.1 


52.3 


2.720 X 10*° 


1.668 X 10'^ 


''^Ca Ar -1- 6n - 26" + 2v^ 


0.42 


8.5 


57.7 


3.145 X 10*° 


2.016 X 10'^ 


62 Ar -^«« s + 6n - 2e- 2v^ 


0.47 


9.0 


63.7 


3.723 X 10™ 


2.488 X 10'^ 


^^'S> ^™ Si + 6n - 2e- + 2v^ 


0.53 


9.4 


70.5 


4.549 X 10™ 


3.153 X 10'^ 


5°Si ^^'^ Mg -1- 8n - 2e- 2v^ 


0.61 


8.8 


79.0 


4.624 X lO''" 


3.472 X lO"' 


"''Mg -^^^ Ne -1- 6n - 2e- -h 2v^ 












*^Ne -^*'^ Ne -^'''^ Ca 


0.66 


10.6 


251.8 


5.584 X lO-"" 


4.399 X 10'^ 


'^Ca Ar + 6n - 26" + 2v^ 


0.69 


4.8 


25.3 


6.883 X 10*° 


5.355 X 10'^ 


66 ^60 s _^ e„ _ 2e- 2v^ 


0.72 


4.7 


27.3 


8.749 X 10*° 


6.655 X 10'^ 


60s Si -h 6n - 2e- -f 2v^ 


0.75 


4.6 


29.2 


1.157 X lO''' 


8.487 X 10'^ 


""Si Mg -1- 8n - 2e- -f 2v^ 












^^Mg -K**^ Mg ^''^ 


0.79 


4.0 


139.6 


1.234 X lO''" 


9.242 X 10'^ 


«^Cr Ti -h 6n - 2e" 2^e 


0.80 


2.0 


8.9 


1.528 X 10*^ 


1.096 X 10'* 


86rp. ^80 ^ g„ _ 2e- + 2!ye 


0.82 


1.9 


9.0 


1.933 X 10*^ 


1.317 X 10'* 


*°Ca -^'^^ Ar -1- 6n - 26" -h 2v^ 


0.83 


1.8 


8.8 


2.510 X 10*^ 


1.609 X lO'* 


■^^Ar S 6n - 2e- -f 2j/e 


0.85 


1.7 


10.2 


3.363 X 10*' 


2.003 X 10'* 


"''S ^''^ Si + 6ri - 2e- -h 2z/e 












62si+62 Si ^'2* Ni 


0.86 


1.7 


70.3 


4.588 X 10*' 


2.520 X 10'* 


"'"Ni -»"'° Fe -h 4n - 26" -h 2i/e 


0.87 


0.8 


2.6 


5.994 X 10*^ 


3.044 X 10'* 


i20pe ^118 2n - 2e- 2:/^ 


0.88 


0.9 


2.4 


8.408 X 10*' 


3.844 X 10'* 


iiscr ^'"5 Ti 2n - 26" -|- 2iye 


0.88 


0.8 


2.2 



the nearly same nuclei after the first pycnonuclcar reac- 
tion, with the slight shift in the densities of the boundaries 
between subsequent shells. The total energy per nucleon 
released then above the neutron drip point is between 1.1 
and 1.5 MeV, with unchanged crustal heating below the 
neutron drip. The value obtained for A\ = 56 is slightly 
larger than that quoted in HZ, which results from correct- 
ing too large neutrino losses in electron captures assumed 
in HZ. 

4. Discussion and conclusion 

We studied nuclear composition of accreted neutron-star 
crusts, assuming different compositions of the X-ray bursts 
ashes. Obtained compositions {A, Z) of the outer crust 
for ashes with = 56 and Ai ~ 100 keep the difference 
by a factor of about two up to the neutron-drip point. 
However, total crustal heating in the outer crust is similar 



for both Ai, and is negligible compared with the crustal 
heating in the inner crust. A chain of processes occurring 
after the neutron drip leads to convergence of composi- 
tions to a common one, which at densities higher than 
10^^ g cm"-^ does not depend on Ai. While the number 
of the heat-sources depends on ^i, and is significantly 
larger for Ai ~ 100 than for Ai = 56, total crustal heat- 
ing, 1.1 MeV/nucleon and 1.5 MeV/nucleon, is quite sim- 
ilar. Simulations with other values of Ai confirm that the 
crustal heating is only moderately dependent on the nu- 
clear composition of the X-ray bursts ashes. 

Let us mention, that in a very recent paper Schatz et 
al. (2003) suggest that rp-processes ashes can burn explo- 
sively at density ~ 10^ g cm~^, giving rise to the X-ray 
superburst. Such a burning will lead to photodisintegra- 
tion of the A ~ 100 nuclei and shift composition of the 
ashes toward iron, so that in this case the original scenario 
of HZ will hold. 
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In view of the importance of the crustal heating for the 
SXRTs models, more detailed study of non-equilibrium 
processes in accreting neutron-star crusts is desirable. An 
investigation along these lines is now being carried out 
and its results will be presented in a separate publication 
(Haensel and Zdunik, in preparation). 
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